Non-axisymmetric Anisotropy of Solar Wind Turbulence 
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A key prediction of turbulence theories is frame-invariance, and in magnetohydrodynamic (MHD) 
turbulence, axisymmetry of fluctuations with respect to the background magnetic field. Paradoxi- 
cally the power in fluctuations in the turbulent solar wind are observed to be ordered with respect 
to the bulk macroscopic flow as well as the background magnetic field. Here, non-axisymmetry 
across the inertial and dissipation ranges is quantified using in-situ observations from Cluster. The 
observed inertial range non-axisymmetry is reproduced by a 'fly through' sampling of a Direct 
Numerical Simulation of MHD turbulence. Furthermore, 'fly through' sampling of a linear superpo- 
sition of transverse waves with axisymmetric fluctuations generates the trend in non-axisymmetry 
with power spectral exponent. The observed non-axisymmetric anisotropy may thus simply arise as 
a sampling effect related to Taylor's hypothesis and is not related to the plasma dynamics itself. 
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Solar wind fluctuations observed by satellites in-situ 
exhibit power law scaling regions identified with an iner- 
tial range of magnetohydrodynamic (MHD) turbulence, 
and with a 'dissipation' range below ion kinetic scales, 
providing a natural laboratory for plasma turbulence (for 
a recent review see, e.g., Ref. In hydrodynamic tur- 
bulence, any anisotropy in fluctuations at large scales 
will tend to isotropize as the cascade proceeds to smaller 
scales 2} . The situation is different in plasma turbulence 
where the existence of a mean magnetic field sets a nat- 
ural preferential direction for anisotropy. Anisotropy is 
thus a key topic in theoretical 0-0], numerical and 
observational studies of plasma turbulence in the solar 
wind [1^3. 

The seminal study of Belcher and Davis [1] used 
Mariner 5 observations to investigate anisotropy of the 
solar wind magnetic fluctuations in the low frequency 
(energy containing) and inertial intervals. They found 
that the fluctuations on average have 5:4:1 power 
anisotropy in an orthogonal coordinate system whose axis 
are [eg xej^g x (eg x e#), eg], where is a unit vec- 
tor in the average magnetic field direction and en is a 
unit vector radially away from the sun. This conclusion, 
that solar wind fluctuations are non-axisymmetrically 
anisotropic with respect to the magnetic field direction in 
the low frequency and inertial intervals was confirmed by 
different authors [3] ■ Recent results using fc-filtering 
were consistent with the results of [8j in that the main 
power was found in the plane perpendicular to the local 
magnetic field distributed preferentially in the direction 
perpendicular to both the magnetic field and the solar 
wind velocity. Dissipation range magnetic fluctuations 
have also been found to be non-axisymmetric [151 ] using 
minimum variance analysis [lfjj. 

The anisotropic expansion of the solar wind can intro- 
duce a preferred direction as captured by models T3j_ 18 1 
and as observed on longer timescales (5-12 firs, see [id]). 



However, from the perspective of turbulence, ordering 
of the observed non-axisymmetric power anisotropy with 
the direction of the solar wind bulk flow velocity at 
the inertial and dissipation scales is rather unexpected. 
If the macroscopic bulk flow speed is sufficiently large 
compared to that of the fluctuations and that of the 
characteristic wave speeds of the plasma, then on the 
timescales over which we observe turbulence in-situ this 
bulk flow simply acts to advect the fluctuating plasma. 
This is Taylor's hypothesis .20;], and if it holds, then 
since the observed properties of the evolving turbulence 
are frame independent they should not correlate with 
the macroscopic flow direction. The observation of non- 
axisymmetry [H, [l(]| with respect to the macroscopic flow 
direction in the high speed solar wind flow is thus para- 
doxical. In plasma turbulence, one would anticipate ax- 
isymmetric anisotropy ordered with respect to the local 
magnetic field. Indeed, theories of MHD and kinetic 
range turbulence assume axisymmetry of statistical char- 
acteristics 21.J221] (note that although the model 
developed in Ref. [4| implies local non-axisymmetry of 
turbulent eddies, it still assumes an axisymmetric en- 
ergy spectrum). As a consequence, studies of anisotropy 
of solar wind turbulence using single spacecraft obser- 
vations often assume axisymmetry (see, e.g., (9L \l2\ l23l| 
and references therein). Understanding the origin of this 
non-axisymmetry is the subject of this Letter, and is es- 
sential if solar wind observations are to be employed in 
the study of turbulence, in particular in the context of 
direct comparisons between theoretically predicted and 
observed statistical properties and scaling exponents. 

Here, we show that the observed non-axisymmetric 
anisotropy can arise as a data sampling effect rather 
than as a physical property of the turbulence. We 
first sample the output of a direct numerical simulation 
(DNS) of MHD turbulence with a 'fly through' emulat- 
ing single spacecraft in-situ observations using Taylor's 
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hypothesis. We will see that this is sufficient to re- 
produce the observed non-axisymmetry in the inertial 
range of the solar wind. To understand how this non- 
axisymmetry can arise, we consider the simplest scenario- 
a 'fly through' sampling of a fluctuating field composed 
of linearly superposed transverse waves with axisymmet- 
ric power anisotropy. The only free parameter in this 
model is the power spectral exponent of perpendicular 
fluctuations. This model reproduces the observed trend 
- that the non-axisymmetry increases with the perpen- 
dicular power spectral exponent as we move from the 
inertial to the dissipation range of scales. The observed 
non-axisymmetric anisotropy may thus simply arise as a 
sampling effect related to Taylor's hypothesis. 

We present the analysis of a sample interval [January 
20, 2007, 1200-1315 UT] of fast quiet solar wind ob- 
served by Cluster spacecraft 4 whilst the magnetic field 
instruments FGM and STAFF-SC were in burst mode, 
providing a simultaneous observation across the inertial 
and dissipation ranges. FGM (sampled at 67 Hz) and 
STAFF-SC data (sampled at 450 Hz) are combined by 
the same procedure as in [2^, [24[ , where a discrete wavelet 
transform is applied to both instrument data sets. This 
merging procedure generates one time series containing 
frequencies ranging from the highest frequency of the 
STAFF-SC data and the lowest frequency of the FGM 
data. This interval is of fast solar wind with a flow speed 
of ~ 590 km/s with plasma parameters: average mag- 
netic field B ~ 4 nT, proton plasma j3 ~ 1.5, proton 
density p p ~ 2 cm -3 , proton temperature T p ~ 29 eV 
and Alfven speed Va — 60 km/s. 

We use the continuous wavelet transform (CWT), as 
outlined in [2|| [26j], to select fluctuations on a specific 
scale, r. The fluctuations are resolved at each r by a 
CWT performed on each component of the magnetic field 
data, B(tj), using the Morlet wavelet, to give a fluctu- 
ation vector, <5B(fj,r). The vector fluctuation are then 
projected onto the local field. At each scale r the local 
magnetic field is defined for every time, tj, by the con- 
volution of a Gaussian window of width 2r centred on 
tj with the data, such that B(tj,r) = [B(tj) * g(tj,r)] 
where g is the Gaussian window. The scale r is related 
to frequency /, (in Hz) of the central frequency of the 
Morlet wavelet. This allows the local magnetic field and 
the fluctuations to be rewritten as a functions of time 
and frequency B(tj, /) and £B(ij, /), respectively. 

We define the local system of unit vectors following Q . 
The unit vector in the direction of the local magnetic 
field is e z (tj,f) = B(tj, /)/ \&(tj, f)\. The other two 
perpendicular unit vectors are ordered with respect to 
the macroscopic flow velocity direction, such that 
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tion time-averaged over the interval. During this interval 
the spacecraft is in a fast and steady stream of the solar 
wind, thus the local velocity is close to the time averaged 
macroscopic velocity in (JT]). 

It can be shown that the results obtained are not sen- 
sitive to the small level of observed variations in the solar 
wind direction. The magnetic field fluctuations are then 
projected to the new basis (P) and the Power Spectral 
Density (PSD) of the corresponding components are de- 
fined by: 



PSD 
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where V is the unit vector of the bulk flow velocity direc- 



where N is the sample size at each frequency and A is 
sampling interval of the data. 

We plot the PSD of these components for this Clus- 
ter interval in Fig. [TJ In the Figure, PSD Z , PSD y and 
PSD X are (from bottom to top) represented by black, 
red and blue lines, respectively. The frequency range 
captured by this interval of data covers both the dissi- 
pation range and the high frequency part of the inertial 
range. The average ratio of the total perpendicular to 
the parallel PSD in the inertial range does not vary sig- 
nificantly with frequency and is ~ 17 : 1. This is con- 
sistent with prior observations [bij that the power in the 
inertial range is predominately perpendicular to the lo- 
cal magnetic field, implying Alfvenic fluctuations. In the 
dissipation range the power isotropises [27], . The ra- 
tio of the PSD of the perpendicular components is shown 
in the inset of the Figure. The mean value of the ratio of 
the PSD of the perpendicular components in the inertial 
range is PSD y /PSD x ~ 0.8, in close agreement with 
This ratio increases significantly in the dissipation range 
to PSD y /PSD x ~ 0.45, coincident with a steepening of 
the PSD power spectral exponents. We will explore the 
origin of these average values and trends. Ratios of the 
PSD can emphasize fluctuations in the PSD that are lo- 
calized in frequency, these can be observational artefacts 
such as spacecraft spin tones, or can be of physical ori- 
gin such as power enhancements at the cross-over from 
the inertial to dissipation range. The PSD ratio can be 
seen to oscillate at the transition between the inertial and 
dissipation range due to these effects. We performed the 
above PSD analysis described on a second Cluster inter- 
val; a fast stream (January 30, 2007, 0000-0100 UT) and 
obtained similar results. 

We now see how this non-axisymmetry can arise. Spa- 
tial snapshots of the three-dimensional velocity and mag- 
netic fields of developed incompressible MHD turbulence 
provide a model for the inertial range fluctuations, these 
were obtained from the DNS described in Q. We will 
consider two simulations: Case I is a globally isotropic 
freely decaying turbulence with resolution 512 3 and Case 
II corresponds to 1024 2 x 256 forced turbulence simula- 
tion with strong background magnetic field Bq , such that 
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b rms /Bo ~ 0.2 Q. A 'fly through' sampling is performed 
on these spatial snapshots by sampling a 'time series' of 
one-dimensional data along different straight line paths 
through the simulation domain at constant velocity to 
mimic single satellite observation of the turbulence using 
Taylor's hypothesis. The direction of the straight line 
path defines the pseudo- macroscopic flow direction V in 
the components defined in ([T]). From each 'fly through' 
we can then calculate the time averaged ratio of the 
power in the two perpendicular components in the same 
manner as above for the solar wind data. These simu- 
lations offer two contrasting, and informative, numerical 
experiments. In Case I the background field direction is 
free to vary and a wide range of angles between the local 
background field and the pseudo-macroscopic flow can 
be realized. For Case II the 'fly through' is restricted to 
the plane perpendicular to the applied background mag- 
netic field, to ensure that the local background magnetic 
field is nearly perpendicular to the pseudo-macroscopic 
flow direction- this will provide us with a simple test case 
in what follows. The results are presented in Figure [2] 
The blue dotted line corresponds to the solar wind data 
shown in the inset of Figure 1. The crosshatched and 
hatched areas indicate the range of values obtained from 
a number of simulation 'fly throughs' for Case I and Case 
II, respectively. Looking at the inertial range, we then 
see a remarkable agreement between the average of the 
power ratio in the inertial range in the solar wind, and 
Case I. Case II also shows non-axisymmetry but tends to 
overemphasize the value as compared to that observed 
by Cluster in the solar wind. 

A full understanding of the origin of power non- 
axisymmetry in the solar wind should capture both that 
it is seen in the inertial and dissipation range of scales 
and that it is strongest in the dissipation range. We now 
show that these features can be qualitatively captured 
by quite simple considerations. Consider a linear super- 
position of waves transverse to a constant background 
magnetic field Bo- For simplicity we assume V SUJ _L Bq. 
To fix coordinates, the solar wind velocity V su , and Bo 
are directed along y and z, respectively so that the mag- 
netic fluctuations <$B associated with the waves are in 
the x, y plane. If the Fourier amplitudes of the fluc- 
tuations are <$B(k) then under Taylor's hypothesis the 
energy densities E y (uj) and E x (oj) observed for the com- 
ponents parallel and perpendicular to V slu are given by 
(see also ^) 

E vM = J d 3 k|<5B, i;c (k)| 2 5( w - k • V sw ). (3) 

If we now assume that there exists some scaling relation 
between the parallel and perpendicular wave numbers 
(for example, critical balance [3() and given 5B(k) _L k, 
after integration over k|| we obtain 

E y (w) = [ d 2 k ± E 2D (k ± )S(u> - k y V sw )sm 2 a, (4) 



E x (u) = j d 2 k ± E 2D (k ± )5{uj - k y V sw )cos 2 a, (5) 

where E 2 u(kj_) is two-dimensional spectrum of fluctua- 
tions [3[ and a is the angle between k^ and V SUJ . Ex- 
pressions ((4]) and ([5]) generally integrate to give E y (ui) ^ 
E x (ui). In particular, for axi-symmetric fluctuations 
E 2 o{k±_) = E 2 D{ki_) = CkJ 1 ^ 1 , where 7 corresponds 
to a one-dimensional spectrum, Eqs. (jU[5]) yield 

T-3 

E y {u) = C^- f(^r + k 2 ) 2 dk x , (6) 

* sw J \ * sw / 

E x {u) = C J k 2 x {^r + k 2 )j 2 dk x . (7) 

This can be integrated in terms of Beta functions and 
after straightforward manipulation (for 7 > 0) we obtain 

E v (u)/E x (u>) = I/7. (8) 

Thus the ratio E y {uj) / E x {uj) is a decreasing function 
of the spectral index in qualitative agreement with 
the results presented above. The values of the ra- 
tio E y (oj)/E x (oj) for different values of 7 are indicated 
on Figure 2. Specifically, for 7 = 1.5 ± 0.1 (iner- 
tial range spectral index observed here) E y (w) / E x (co) = 
0.67 ± 0.04 (gray area in the low frequency range), for 
7 = 5/3 (Kolmogorov spectrum) E y (ui) / E x (ui) = 0.6 
(black solid line), for 7 = 3/2 (Iroshnilov-Kraichnan 
spectrum) E v (uj) / E x (uj) — 0.67 (red dotted line) and for 
7 = 2.76±0.05 (dissipation range spectral index observed 
here) E y (uj) / E x (u>) — 0.36 ± 0.02 (gray area in the high 
frequency range). 

This result has a simple physical explanation. When 
the satellite samples three-dimensional fluctuations as 
a one-dimensional series all the fluctuations with fixed 
k y contribute to a given frequency w = k y V sw . This 
mixture of different modes does not produce any non- 
axisymmetric anisotropy if energy density is equally dis- 
tributed among different scales (7 = 1), but if this is 
not the case, the resultant frequency spectra for the 
energy density of different components will be non- 
axisymmetric. Indeed, consider the case when energy 
density is decreasing with scale (7 > 1). Then, at a 
given frequency / = k y V sw , since the waves are trans- 
verse, modes with k x < k y contribute more to E x (uj) 
compared to E y (u>). For the modes with k x > k y the 
situation is reversed, but because the modes in the latter 
range have less power, one finally has E x (u>) > E y {ui). 

Results from the simple model can be seen from Fig- 
ure 2 to be consistent with the results from Case II of 
the DNS which share the same restricted geometry, that 
is, the background field is perpendicular to the macro- 
scopic flow velocity. The model also gives the qualitative 
trend in the anisotropy with the spectral index of the tur- 
bulence. On the other hand, the simple model predicts 
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stronger power anisotropy than seen in the solar wind as 
well as in the DNS without strong background magnetic 
field (Case I). This suggests that the overestimation in 
non-axisymmetry is indeed due to this restricted geom- 
etry of the simple model. To practically test this idea 
would require knowledge of the distribution of the angle 
between V su , and B which, within any given interval of 
solar wind data is non-uniform, and which can vary from 
one interval to another. 

The sampling effects described above have implications 
for the results obtained by methods used for the study 
of turbulence anisotropy that rely upon Taylor's hypoth- 
esis. In addition, multispacecraft methods that do not 
use Taylor's hypothesis explicitly, such as fc-filtering (see, 
e.g., [10( and references therein) still require filtering of 
the data in the spacecraft frame - selecting perturbations 
within some frequency interval [— u> ma x , w mm ] . As the so- 
lar wind speed is much larger than the Alfven speed, it is 
clear that this procedure is essentially a filtering with re- 
spect to the wave vector component parallel to the solar 
wind velocity. As demonstrated here, this leads to pref- 
erential filtering of the modes with wave vectors mainly 
parallel to the solar wind, and as a result could generate 
non-axisymmetric anisotropy. 

Finally, we have for the restricted case where the back- 
ground field is perpendicular to the flow, obtained an 
explicit relationship between the spectral exponent and 
the power non-axisymmetry. For this special case, the 
non-axisymmetry in power is quite sensitive to the value 
of the exponent. It may be possible to develop this for- 
malism, at least numerically, to use the observed non- 
axisymmetry in power to determine the scaling exponent 
more accurately. This is central to testing predictions of 
turbulence theories. 
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Figure 1: PSDs of magnetic field components (from bottom 
to top) - PSD Z (black), PSD y (red) and PSD X (blue), where 
e* = B/ bI , e x = e z x V/ |e z x V| and e y — e z x e^. Error 
bars are alwavs smaller than 4% and are usually smaller than 




Figure 2: The ratio E y /E x , where and refer to a field 
aligned coordinate system as in Figure [T] for the solar wind 
data (blue dotted line), the range of values obtained from a 
number of simulation 'fly throughs' for Case I (crosshatched 
area) and Case II (hatched area), as well as predictions of 
Eq. ||HJ| for the spectral indices 7 = 5/3 (black solid line), 
7 = 1.5 (red dashed line), 7 = 1.5 ± 0.1 (left gray area) and 
7 = 2.76 ± 0.05 (right gray area). 



